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Pion transverse charge density from timelike form factor data 



O 

> 
O 

in 

Oh 1 

<D 

43 



> 
(N 



X 



G. A. Miller, 1 M. Strikman, 2 and C. Weiss 3 

'University of Washington, Seattle, WA 98195-1560, USA 
2 Department of Physics, Pennsylvania State University, University Park, PA 16802, USA 
Theory Center, Jefferson Lab, Newport News, VA 23606, 
(Dated: November 2, 2010) 



USA 



The transverse charge density in the pion can be represented as a dispersion integral of the imag- 
inary part of the pion form factor in the timelike region. This formulation incorporates information 
from e + e~ annihilation experiments and allows one to reconstruct the transverse density much 
more accurately than from the spacelike pion form factor data alone. We calculate the transverse 
density using an empirical parametrization of the timelike pion form factor and estimate that it is 
determined to an accuracy of ~ 10% at a distance h ~ 0.1 fm, and significantly better at larger 
distances. The density is found to be close to that obtained from a zero-width p meson pole over a 
wide range and shows a pronounced rise at small distances. The resulting two-dimensional image of 
the fast-moving pion can be interpreted in terms of its partonic structure in QCD. We argue that 
the singular behavior of the charge density at the center requires a substantial presence of pointlikc 
configurations in the pion's partonic wave function, which can be probed in other high-momentum 
transfer processes. 

PACS numbers: 11.55.Fv, 13.40.Gp, 13.60.Hb, 13. 66. Be 
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I. INTRODUCTION 

Learning to describe the structure and interaction of 
hadrons on the basis of QCD is one of the main objectives 
of nuclear physics. An essential step in this program is 
to understand the structure of the pion, a nearly mass- 
less excitation of the QCD vacuum with pseudoscalar 
quantum numbers. The pion plays a central role in nu- 
clear physics as the carrier of the long-range force be- 
tween nucleons and a harbinger of spontaneous symme- 
try breaking. The importance of the pion has been rec- 
ognized by intense experimental and theoretical activity 
aimed at measuring its properties and understanding its 
structure. The pion electromagnetic form factor F K (t) 
was measured at spacelike momentum transfers through 
pion-electron scattering 0, 0| and pion electroproduc- 
tion on the nucleon [3|-|6|; new measurements in the re- 
gion \t\ ~ few GeV 2 are planned with the Jefferson Lab 
12 GeV Upgrade 0. In the timelike region the modu- 
lus of the (complex) pion form factor, lfir(£)L was deter- 
mined in a series of e + e~ experiments [1413; see R e f- 03 
for a compilation of the older data. 

The concept of transverse densities 0], whose prop- 
erties were explored in several recent works (lH [H| , pro- 
vides a model-independent way to relate the form factors 
of hadrons to their fundamental quark/gluon structure in 
QCD. Defined as the 2-dimensional Fourier transforms of 
the elastic form factors, the transverse densities describe 
the distribution of charge and magnetization in the plane 
transverse to the direction of motion of a fast hadron; see 
Ref. [I?) for a review. They are closely related to the par- 
ton picture of hadron structure in high-energy processes 
and correspond to a reduction of the generalized par- 
ton distributions (or GPDs) describing the distribution 
of quarks/ antiquarks with respect to longitudinal mo- 



mentum and transverse position [l8l [l9j . It is therefore 
natural to attempt to interpret the pion form factor data 
in terms of the transverse charge density in the pion. 
In particular, the density at small transverse distances 
b <C 1 f m places constraints on the probability of point- 
like configurations (or PLCs) in the pion — qq configu- 
rations in the partonic wave function of a transverse size 
much smaller than the typical hadronic radius [20] . Such 
configurations play an important role in high-momentum 
transfer reactions involving pions, such as the pion tran- 
sition form factor 7*7 — > n° I2ll[22l o r pion production in 
large-angle scattering processes |23{. They are essential 
for the physics of the color transparency phenomenon 
predicted by QCD [H, [25|, which is studied in high- 
energy pion dissociation on nuclear targets [H, [2?} and 
electromagnetic pion knockout [H, [2i| and is closely re- 
lated to the existence of factorization theorems for hard 
meson production processes. The dynamical origin of 
PLCs — whether they are generated through perturba- 
tive QCD interactions with large-size configurations or 
by non-perturbative mechanisms, remains a subject of 
intense study. 

The transverse charge density in the pion is defined 
as the 2-dimensional Fourier transform of the spacelike 
pion form factor, 
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where F v is regarded as a function of the invariant mo- 
mentum transfer t. The function p^(b) gives the prob- 
ability that charge is located at a transverse separa- 
tion b from the transverse center of momentum, with 
f d 2 b pir{b) — 1. The definition Eq. ([T]) may in princi- 
ple be used to calculate the charge density directly from 



2 



the spacelike form factor data. In the nucleon case, 
where the spacelike form factors can be extracted di- 
rectly from the measured eN elastic scattering cross sec- 
tion and are known up to rather large momentum trans- 
fers, this approach has been quite successful; see Ref. [30( 
for an assessment of the uncertainties. In the pion case 
the spacelike form factor at momentum transfers above 
Q 2 > 0.25 GeV 2 was extracted only indirectly in electro- 
production experiments on the nucleon iV(e, e'w)N' , with 
substantial model dependence, and is known only poorly 
at higher Q 2 , rendering such a program difficult. How- 
ever, for the pion one has another avenue for evaluating 
the transverse density, based on a dispersion representa- 
tion for the pion form factor. Noting that the singulari- 
ties of F„ (t) as an analytic function of t are confined to a 
cut along the positive real axis starting at t — 4m 2 , the 
form factor can be expressed as [31| 



dt' 



Im^(i') 



t' -t + iO 
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The asymptotic behavior expected from perturbative 
QCD, F n (t) ~ a s (t)/\t\ for t — >• oo, allows the use of 
an unsubtracted dispersion relation 14411. Substitution of 
Eq. J2J in Eq. (P leads to the result [32] 
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This representation of the charge density as a dispersion 
integral over the imaginary part (or spectral function) of 
the timelike pion form factor has an interesting "filter- 
ing" property. The exponential drop-off of the modified 
Bessel function Kq at large arguments causes the inte- 
grand of Eq. to decrease exponentially at large t and 
ensures that only values y/i ~ 1/b in the spectral func- 
tion are effectively sampled at a given distance b. In 
the nucleon case the timelike form factor is measurable 
only at t > Am 2 N and Eq. ([3]) is not useful for calcu- 
lating the transverse density from data (it is, however, 
very useful for theoretical analysis; for example, the chi- 
ral large-distance component of the nucleon charge den- 
sity at b ~ m^ 1 can be obtained from the calculable 
strength of the two-pion cut in the nucleon form factor 
near threshold [32]). In the pion case the physical region 
for the timelike form factor starts at t = 4 m 2 , covering 
the entire range of the dispersion integral, and Eq. (|3]) 
becomes a practical method for calculating the charge 
density at all values of b. High-quality e + e~ annihila- 
tion data exist for values of t up to ~ 1 GeV 2 , so that we 
hope to be able to determine p-n(b) accurately for values 
of b at least down to values of b ~ 1 GeV -1 = 0.2 fm. 

The imaginary part of the pion form factor 1m.F K (t) 
entering in the dispersion representation Eq. is not 
measured directly in annihilation experiments. The 
e + e _ — > 7r + 7r~ cross section is proportional to \Fn{t)\ 2 , 



and model-dependent input is generally needed to de- 
termine the phase. In the region of the p meson reso- 
nance this problem was studied extensively long ago and 
is under good theoretical control. The phase of the first 
higher resonance p' is strongly constrained by the dis- 
persion integrals (sum rules) for the pion charge and the 
measured charge radius. At larger values of t arguments 
based on perturbative QCD and local duality provide 
some guidance. Combined with the filtering property of 
the dispersion integral Eq. , these constraints strongly 
reduce the model dependence in the transverse density at 
6 > 0.1 fm. Our estimates below show that the this way 
of constructing p v (b) gives substantially more accurate 
results than use of the spacelike pion form factor data 
alone. 

In this article we calculate the transverse charge den- 
sity in the pion in the dispersion representation Eq. ([3]) 
using an empirical parametrization of the timelike pion 
form factor based on e + e~ annihilation and spacelike 
form factor data [33| . We find that the density is deter- 
mined to an accuracy of ~ 10% at transverse distances 
b ~ 0.1 fm, and substantially better at larger values. We 
thus obtain a precise two-dimensional image of the fast- 
moving pion, which can be interpreted in terms of its 
partonic structure in QCD. In particular, the density 
exhibits a pronounced rise at small &, as was observed 
earlier — although with much lower precision — in an 
analysis based on the spacelike pion form factor [l6j . Us- 
ing experimental information on the quark density in the 
pion, we argue that such singular behavior of the charge 
density cannot be explained by large-size, x — > 1 config- 
urations in the pion's partonic wave function and must 
therefore be attributed to PLCs. Our result thus places 
constraints on the probability of PLCs in the pion, which 
can be probed in other high momentum-transfer pro- 
cesses involving pions. 

The plan of this paper is as follows. In Sec.|H]we briefly 
describe the main features of the pion form factor in 
the timelike region and the elements of the parametriza- 
tion of Ref. [3j]. In Sec. IIIII we calculate the transverse 
charge density and investigate its uncertainties at small 
distances. The implications for the pion's partonic struc- 
ture and the presence of PLCs are discussed in Sec. IIVI 
Section [V] discusses the possible role of chiral dynamics 
in the pion transverse density at large distances. A sum- 
mary and suggestions for further studies are presented in 
Sec. ED 



II. TIMELIKE FORM FACTOR 
PARAMETRIZATION 

In the energy region y/i < 1 GeV the measured pion 
form factor |F w (t)| 2 is dominated by the p meson res- 
onance, with clearly visible effects of p-cu mixing (see 
Ref. [33| for a summary of the data). Theoretical sup- 
port for p dominance at the amplitude level comes from 
the observation that the 2tt channel accounts for most of 
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the annihilation cross section, which allows one to relate 
the pion form factor to the 7T7t scattering amplitude via 
elastic unitarity. In this region the form factor is success- 
fully described by the Gounaris-Sakurai (GS) amplitude 
34], which is derived from an effective range expansion 
of the 7T7T phase shift and has the correct analytic struc- 



ture. The neglect of certain off-shell terms oc (t 
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in the GS amplitude leads to a Breit-Wigner (BW) type 
parametrization with energy-dependent width; this sim- 
plified form also describes the |i^(t)| 2 data in the region 
\/t ;$ 1 GeV but does not respect the analytic properties 
of the form factor (it has a spurious branch cut singularity 
at t — 0). We shall employ the full GS parametrization 
in our studies here. 

Above the p region data for |F,,.(i)| 2 exist up to ener- 
gies \ft < 3 GeV. Because of the many hadronic channels 
in the total cross section, the phase of the form factor at 
these energies is much more uncertain. In the region 
of the first higher resonance p' the phase is constrained 
by the sum rules for the pion charge and the charge ra- 
dius, which require partial compensation of the spectral 
strength in the p meson region. At higher energies the- 
oretical constraints come from the asymptotic behavior 
predicted by perturbative QCD, which demands strong 
cancellations between higher resonances in a resonance- 
based description, as indeed found in dual resonance 
models. 

In the present study we use the timelike pion form 
factor parametrization of Ref. [33[, which describes the 
high-energy region by a pattern of resonances consis- 
tent with the QCD asymptotic behavior. The param- 
eters were determined by a detailed analysis of the time- 
like data up to yt < 3 GeV. The continuation of these 
parametrizations to t < also describes the spacelike 
form factor in accordance with the data, including the 
recent JLab data up to \t\ — 2.45 GeV 2 p|, which ap- 
peared after publication of Ref. [33j . 

A brief description of the elements of the parametriza- 
tion of Ref. [33| is provided here; for details we refer to 
the original article and references therein. The first four 
p meson resonances are included as specific states with 
masses up to 2.0 GeV (p-u mixing is taken into account 
for the lowest resonance) . These resonances are described 
by the GS form, which incorporates the proper threshold 
behavior of the widths and has the correct analytic prop- 
erties [45l | . In addition, an infinite series of higher excita- 
tions is included via an ansatz [351 ] based on the dual reso- 
nance model. Its continuation to the space-like region ex- 
hibits a smooth behavior with a power-law asymptotics 
as with (3 — 2.1 — 2.3. The imaginary part of the 

form factor obtained with the GS parametrization [33j is 
shown in Fig. QJi (solid line) . One clearly sees the domi- 
nance of the p meson pole in the region \ft < 1 GeV, and 
the alternating sign of successive resonance contributions 
at larger values of y/t, as expected from theoretical con- 
siderations. 

To estimate the uncertainty in the imaginary part, we 
have taken the quoted variances of the fit parameters of 
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FIG. 1: (a) Solid line: The imaginary part of the pion form 
factor obtained from the fit of Ref. [33] (GS parametrization), 
as a function of i/i. Shown here is the function \fi ImF n (t)/ir, 
which effectively enters in the dispersion integral over yi, 
Eq. ((3|. Dotted lines: ±ler error resulting from the uncorre- 
cted uncertainties of the fit parameters. The threshold en- 
ergy yt = 2m,r is indicated by the vertical line, (b) The 
weight factor Ko(Vtb) in the dispersion representation of the 
transverse charge density Eq. (|3} , as a function of yt, for sev- 
eral values of b. Shown are the functions normalized to unity 
at the threshold \ft — 2m^. 



Ref. [33| and studied the statistical variation of the imag- 
inary part, assuming uncorrelated errors. The resulting 
±1<7 error band is shown in Fig. [TJi (dotted lines). The 
variance in the p meson mass region is at the few per- 
cent level. At energies above 1 GeV it becomes substan- 
tially larger, reaching close to 100% at yt = 2 GeV. Note 
that in this energy region our uncorrelated estimate likely 
represents an upper bound on the uncertainty, as correla- 
tions between the statistical fluctuations of the coupling 
and width of the second resonance would considerably re- 
duce the overall fluctuations of the imaginary part near 
\ft ~ 1.4 GeV. For energies above 3 GeV we cannot re- 
liably estimate the relative uncertainty of the imaginary 
part in this way, as the couplings of the resonances in this 
region are dictated by the dual resonance model, and the 
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data in this region are very poor. However, the imagi- 
nary part in this region is expected to be very small and 
contributes negligibly to the charge density at b > 0.1 fm 
(see below) , so that its relative uncertainty is not impor- 
tant for our purposes. We emphasize that we use the 
parametrization of Ref. [33] only as an effective repre- 
sentation of Im F v (t) in the energy range y/i < 3 GeV, 
and that our conclusions do not depend on the particu- 
lar y/i — > oo asymptotic behavior imposed by the dual 
resonance model. 



III. TRANSVERSE DENSITY AND ITS 
UNCERTAINTY 

We now use the timelike form factor parametrization 
to evaluate the transverse charge density in the pion and 
estimate its uncertainty. It is instructive to study first the 
distribution of strength in the dispersion integral Eq. ([3]) . 
The imaginary part Im F w (t) is weighted with the mod- 
ified Bessel function K^y/ib), which exponentially sup- 
presses energies y/i ^> 1/6. Fig. [Tb shows this weight 
factor as a function of y/i for several values of b, normal- 
ized to the same value at threshold y/i = 1m^ 1 i.e., the 
ratio 
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K (Vib)/Ko(2m„b). 



(4) 



One sees that the effective distribution of strength in y/i 
strongly changes with the distance b. At b = 0.1 fm a 
noticeable contribution to the dispersion integral comes 
from the region y/i > 1 GeV, where the parametrization 
of ImF,(t) shows considerable uncertainty (see Fig.[T^). 
At b = 0.5 fm these contributions are largely suppressed, 
resulting in almost perfect "vector meson dominance" in 
the dispersion integral. Finally, going to distances as 
large as b ~ 2 fm, one begins to suppress also the p mass 
region and emphasizes the near-threshold region of the 
form factor, y/i — 2771^ ~ fewm,r. 

In order to quantify the accuracy of the calculated 
transverse density we need to study the numerical con- 
vergence of the dispersion integral at large values of y/i. 
Figure [5] shows the percentage deviation of p^ (b) from 
the full result as a function of a cutoff applied to the up- 
per limit of the y/i integral in Eq. ((3]) (here the integral is 
evaluated with the central value of the GS parametriza- 
tion as shown in Fig. [T]). One sees that at b — 0.1 fm 
the region y/i > 3 GeV accounts for only about ~ 1 % of 
the total integral, meaning that even a drastic change of 
Im F„(t) in this region by a factor 2 — 3 would change the 
density only by ~ 2 — 3% (46|. The error in the density 
is thus dominated by the mass region 1 < y/i < 3 GeV, 
where we have estimated the uncertainty of Im F v (t) in 
Sec. HI With a - 100% uncertainty at y/i = 2 GeV, 
where the integral has converged to within ~ 4% of its 
value, we expect an uncertainty of the density of (at least) 
~ 4%. Surprisingly, even for much smaller distances the 
region y/i > 3 GeV seems to contribute relatively little 
to the dispersion integral; see the curve in Fig. [5] for 
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FIG. 2: Percentage deviation from the full result for the dis- 
persion integral Eq. @, as a function of the upper limit 
of y/i, for b — 0.5 fm (solid line), 0.1 fm (dashed line) and 
0.02 fm (dotted line). The integrand is evaluated using the 
GS parametrization of Ref. • 



b ~ 0.02 fm. While the integral requires larger values of 
y/i to converge, the contribution from y/i > 3 GeV is still 
only ~ 2%, and the overall uncertainty can be estimated 
from that of the 1 < y/t < 3 GeV. At larger distances 
b ~ 0.5 fm, however, the integral has fully converged al- 
ready at y/i ~ 1 GeV, and the overall uncertainty is dom- 
inated by the low-energy region y/i < 1 GeV. In this 
region the parameter errors in the fit are so small that 
the model dependence of the parametrizations (details of 
p—oj mixing, p line shape) can no longer be neglected in 
establishing the overall error. 

Given the dominance of energies y/i < 3 GeV in the 
dispersion integral, we can evaluate the density with the 
parametrization of Ref. [33[ and estimate its uncertainty 
from the parameter error band shown in Fig. [T] The result 
is displayed in Fig. [3] The quoted la error in ImF^t) 
translates into an uncertainty of p^ib) of ±(1.5, 7, 13)% 
at b = (0.5, 0.1, 0.02) fm. The density is thus determined 
much more accurately, and down to much smaller dis- 
tances, than from the spacelike pion form factor data 
alone |l6[. 

A welcome feature of the dispersion representation of 
the charge density, Eq. @, is that the kernel K${y/ib) is 
a positive function. As a result, an upper or lower bound 
on the spectral function lmF w (t) directly provide a cor- 
responding bound on p„-(6), greatly simplifying the error 
analysis. (A method to estimate the uncertainty of the 
charge density as the Fourier transform of the spacelike 
form factor was described in Ref. |30(.) 

An additional source of uncertainty in the charge den- 
sity at small distances are recent data on the timelike 
pion form factor at large values of t that were not in- 
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FIG. 3: Transverse charge density in the pion, p v {b). Solid 
line: Dispersion integral Eq.([3} evaluated with the GS form 
factor parametrization [331 ] (see Fig.QJ,). Dashed lines: la er- 
ror resulting from the quoted uncertainty of the parametriza- 
tion (see Fig.QJ,). Dotted line: Density resulting from a heavy 
resonance with mass tur = 3.67 GeV and width Yr = 0.2 thr, 
providing a rough assessment of the impact of the CLEO time- 
like form factor data [1 21 ] (details see text). 

eluded in the fit of Ref. [33j]. The CLEO measurement 
12] at Vt = 3.67 GeV reports a value of \F„\ = 0.075 ± 
0.008 (stat) ± 0.005(syst), much larger than the value 
0.034 provided by the GS parametrization of Ref. [33| . 
We see no simple way to modify the parametrization to 
account for this datum. Indeed, Ref. [33| argues that 
increasing the absolute value of the form factor by a fac- 
tor of ~ 2 at large \ft is not possible. In particular, the 
article states that it is implausible for the form factor ob- 
tained on the basis of a dual resonance parametrization 
to reach values ^(t)^ > 0.01 at y/i = 2.5 - 3 GeV (as 
would correspond to the new datum, assuming power- 
like i-dependence) without conflicting with the spacelike 
data and especially with QCD predictions [47]]. 

One possibility is that the error of the CLEO result is 
larger than estimated in Ref. [13]. Another possibility is 
that there is a new mechanism providing a strong cou- 
pling to two pions at high energies. Here we only wish 
to make a rough assessment of the potential impact of 
this new datum on the transverse density. To this end, 
let us assume the existence of an "additional" irir res- 
onance at \ft = m R = 3.67 GeV, described by the GS 
form, whose coupling c R to the virtual photon is related 
to the measured pion form factor as 

\F w (m R )\ = c R m R /T R . (5) 

Taking the width T R to be ~ 20% of the mass, as it 
is for the p meson, we obtain a coupling c R — 0.015 
from the CLEO measurement. Such an addition gives 
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FIG. 4: Transverse charge density p^ib) in the pion. Solid 
line: Dispersion integral Eq. ([3]) evaluated with the GS form 
factor parametrization of Ref. [331 - Dashed line: Density from 
zero-width p meson pole, Eq. ©. 

a negligible contribution to the 1^(^)1 at values of \fi 
for which most of the data entering in the parametriza- 
tion [H[ were taken; for example, it provides a ~ 1% 
contribution to at -Jt = 1 GeV. The "extra" con- 

tribution to the charge density from such a resonance 
would be +(0.04, 4, 16)% at b= (0.5, 0.1, 0.02) fm (see 
Fig. [3]). If we added this uncertainty to the one esti- 
mated previously from the error of the parametrization 
for \ft < 1 GeV, we would conclude that the density is 
determined to (±1.5, +11 - 7, +39-16) % at the quoted 
values of b. This is surely a conservative estimate, as at 
least part of the uncertainty in the unmeasured high-i 
region is already included in the parametrization error. 
A larger value of the width of the hypothetical resonance 
would lead to a proportionately larger contribution to 
p^{b), but would have to be reconciled with the precise 
data for |F ff | in the mass region \Ji ^ 1 GeV. We con- 
clude that the new CLEO data have only a modest im- 
pact on the transverse density at distances ~ 0.1 fm, but 
may cause substantial modifications at smaller distances. 

Figure [3] shows the transverse density obtained from 
the dispersion integral on a logarithmic scale, which al- 
lows one to see the exponential fall-off at larger distances. 
For comparison we also show the density obtained from 
a single resonance of zero width at the p meson mass m p , 
with a coupling chosen to ensure unit charge (i.e., the 
vector meson dominance model) 

P;r(fr)zoro-width = {m 2 p / 2ir) K (m p b) . (6) 

One sees that the dispersion result is very close to the 
zero-width p form for all distances 0.1 < b < 1 fm 
and can be represented by the latter within the esti- 
mated errors (at larger values of b the spectral strength 
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FIG. 5: (Color online) Three-dimensional rendering of the 
transverse charge density in the pion, as obtained from the 
dispersion integral Eq.lO evaluated with the GS form factor 
parametrization of Ref . [33| : cf. Figs. [3] and [4] 
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FIG. 6: Probability accumulation Eq. ((7]) in the transverse 
density (cf. Figs. [3] and f3]). Solid line: Dispersion integral 
(GS parametrization). Dashed line: Zero-width p meson pole. 
The arrow indicates the experimental RMS transverse charge 
radius. 



near threshold becomes important; see Sec. [V}. What 
is more, the dispersion result follows the zero-width p 
curve down to much smaller distances, being only a few 
percent smaller down to b = 0.01 fm. This shows that 
there are very strong cancellations between the effective 
poles parametrizing the high-mass continuum. As we 
just demonstrated, there is considerable uncertainty in 
the dispersion result for the density at such small dis- 
tances. However, there is the intriguing possibility that 
the density might effectively be described by vector me- 
son dominance down to distances significantly smaller 
than the inverse p meson mass, m^ 1 = 0.25 fm. 

In Fig. [5] we show a 3-dimensional rendering of the 
transverse charge density, which conveys also the infor- 
mation on the supporting area and thus gives an impres- 
sion of the true physical shape of the fast-moving pion 
as seen by an electromagnetic probe. Our dispersion ap- 
proach provides a data-based image of the pion's trans- 
verse structure at small distances with unprecedented 
precision. One clearly sees the strong rise of the trans- 
verse density toward the center. This remarkable obser- 
vation calls for a microscopic explanation in terms of the 
pion's partonic structure. 



IV. IMPLICATIONS FOR PION PARTONIC 
STRUCTURE 

The results of our empirical study of the transverse 
charge density have interesting implications for the par- 



tonic structure of the pion in QCD. The transverse charge 
density puts constraints on the possible distribution of 
transverse sizes of configurations in the pion's partonic 
wave function. A useful quantity to consider is the inte- 
gral of the transverse charge density up to a given dis- 
tance, 

P(b) = j 'd 2 bQ{b~b') p n (b% (7) 

which determines the cumulative probability for configu- 
rations contributing to the transverse density at the dis- 
tance b. The probability obtained from our dispersion 
result for the charge density (cf. Figs. [3] and 2]) is shown 
in Fig. |51 together with that obtained from a zero-width 
p meson pole (cf. Eq. [5]), 

P(&)zero-width = m p bKi(m p b). (8) 

The probability reaches 1/2 at b = 0.33 fm, a value some- 
what smaller than the root of the mean squared (RMS) 
transverse radius, (b 2 )]/ 2 = 0.53 fm. This is to be ex- 
pected, as large-size configurations are counted with a 
higher weight in the average of b 2 than than the me- 
dian. The RMS transverse radius calculated from our 
dispersion integral for the charge density agrees very well 
with the value extracted from the slope of the low-t 
pion form factor measured in ire scattering experiments, 
(r 2 )^ = (3/2)(b 2 ) 7T = 0.439 ± 0.008 fm 2 [J13, as was al- 
ready noted in the discussion of the fit to the timelike 
form factor data in Ref. [33j. 

To understand how the transverse charge density is re- 
lated to the partonic structure it is necessary to recall 
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FIG. 7: Transverse distances in a qq configuration of the 
pion's partonic wave function: b is the distance between the 
quark and the transverse center of momentum, r the distance 
between the q and q. 



the relationship between the coordinate b and the phys- 
ical transverse size of configurations in the fast-moving 
pion. The coordinate b measures the distance between a 
constituent — say, a quark q — and the transverse cen- 
ter of mass of the pion. If the quark carries longitudinal 
momentum fraction x, and the remnant system R carries 
I — x, the transverse center of momentum of the pion is 
at xr q + (1 — x)rn, where r qi R denotes the transverse 
position of the quark and the center of momentum of the 
remnant system. The transverse separation of the quark 
from the remnant system is thus given by 



b/(l-x). 



(9) 



Figure [7] illustrates this relation for a qq configuration 
in which the remnant system consists of a single anti- 
quark. In the transverse charge density one considers 
the charge-weighted density of constituents at a given 6, 
which is obtained as the average over configurations with 
different x and physical size r in the partonic wave func- 
tion. Equation (|9|) now implies that the charge density at 
b much smaller than the typical hadronic size i?h a d ~ 1 fin 
can arise from two different classes of configurations: 

I) Small physical size r <C i?had and non-exceptional 
values of x, i.e., not close to 1 (PLCs). One ex- 
pects the elementary qq configuration accounts for 
a significant fraction of these configurations, as 
the emission of gluons and creation of additional 
qq pairs are suppressed in small-size color-singlet 
configurations. The existence of such PLCs is re- 
quired by the non-zero value of the pion weak decay 
matrix element, parametrized by the constant f^, 
where the axial current operator annihilates a qq 
pair in a point in space [2l| . 

II) Large physical size r ~ i?had and extreme momen- 
tum fractions x ~ f — 6/i?had (end-point configu- 
rations). These are generally not just qq configura- 
tions, as soft gluon radiation is not suppressed in 
large-size configurations. The probability of such 
configurations determines the behavior of the par- 
ton densities in the pion at large values of x. 

As all analysis of partonic structure, the distinction be- 
tween the two classes of configurations depends on the 



resolution scale Q 2 . Standard leading- log Q 2 evolution 
degrades the parton momentum fractions and reduces the 
probability of end-point configurations. The total charge 
density resulting from the sum of all configurations is of 
course scale-independent, being the matrix element of a 
conserved current. 

We can estimate the possible contribution of large- 
size x 1 configurations to the transverse density at 
small b in a simple phenomenological model, using infor- 
mation on the quark distribution in the pion at large x 
extracted from fits to irN Drell-Yan data [36| . Our basic 
assumption here is that the physical transverse size of 
large-x configurations in the pion tends to a finite value 
of the order of the typical hadronic size. Generalizing 
the expression obtained from the overlap of light-cone 
wave functions of individual configurations, we model 
the x— and 6-dependent charge density (i.e., the charge- 
weighted quark GPD) arising from large-a; configurations 
as 



argc — size 



f(r = b/(l-x)) 
q ^ X ' (1-x) 2 ' ( ' 



where q^ (x) is the valence quark distribution in the pion 
48] and f(r) describes the distribution over physical 
transverse sizes r, with a range of the order of the typ- 
ical hadronic size, normalized such that J d 2 r f(r) = 1; 
Eq. (|TU1) thus satisfies J d 2 bp 7r (x,b) = q n (x). The trans- 
verse charge density p v (b) arising from large-size con- 
figurations is then given by the integral of the density 
Eq. (jlOp over x. In calculating this integral we impose 
the physical requirement that the transverse size r of the 
configuration be larger than some critical ro ■ This limits 
the range of x in the integral to values x > l — b/rg, where 
it is assumed that b < r$. We thus consider the "con- 
ditional" large-size contribution to the density defined 
as 



p v {b\r > r ) = 



l-b/r 



dx P-jy (x, ^)largc— size ■ (11) 



To evaluate this contribution to the charge density at 
small 6, we use the the parametrization of the pion quark 
density of Ref. [37j ■ The size distribution f(r) we take to 
be of Gaussian form, f(r) = exp(— r 2 / R 2 )/(ttR 2 ), where 
the parameter R 2 — (r 2 ) n defines the average squared 
radius and is of the order of the typical hadronic size 
~ lfm 2 . For a loosely bound qq state with (x)^ = 1/2 
one would have (r 2 )^ = (b 2 /(l — x) 2 )^ w 4(6 2 ) 7r ; a nat- 
ural choice is therefore R 2 = 4(6 2 ) 7r oxp = 1.16 fin . Fig- 
ure [5] shows the contribution to the charge density from 
configurations with r > ro = 0.2 fin estimated with this 
model, for two values of Q 2 . One sees that it accounts 
only for at most ~ 20% of the total transverse density at 
b = 0.1 fm, and even less at smaller distances. We thus 
conclude that large-size configurations with x — > 1 play 
only a minor role in the pion transverse charge density at 
small b, and that most of it can be attributed to PLCs. 

The small-6 behavior of the large-size contribution to 
the charge density in our model can formally be related to 
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FIG. 8: Dashed/dotted line: Contribution of large-size con- 
figurations with r > ro to the transverse charge density in the 
pion, as estimated in the model defined by Eqs. (|10[) - (|11[1 , for 
ro = 0.2 and two values of Q 2 . Solid line: Density obtained 
from the dispersion integral, cf. Fig. 



the power behavior of the quark distribution in the pion 
for x — > 1. In the limit of small b the integral in Eq. (fTTj) 
extends over a narrow range of x close to 1. If the quark 
distribution vanishes as q- K {x) <~~> (1— x) 13 , one easily shows 
that the density Eq. (fTTj) scales as p n (b\r > r ) ~ fr 3-1 
for b — > 0. The change in the small-6 behavior with Q 2 
seen in Fig. [5] reflects the effect of QCD evolution on 
the exponent j3. Note, however, that even in the low- 
Q 2 region where /? < 1 the large-size contribution in 
our model is substantially smaller than the total density 
obtained from the dispersion integral. 

In sum, our estimate shows that large-size x — ► 1 con- 
figurations cannot account for the strong rise of the trans- 
verse density at small b, and that it is therefore reason- 
able to interpret the empirical density in terms of PLCs 
in the pion's partonic wave function. In qq configura- 
tions of small size, it is expected that the wave function 
peaks at x = 1/2, which implies that the physical trans- 
verse size of the most likely configurations is r 2b. 
With the plausible assumption that the small-size con- 
figurations in the pion are mostly qq, we would conclude 
from Pig. E that there is a probability of 12% (29%) for 
configurations with b < 0.1 fm (0.2 fm), and thus with 
qq separation r < 0.2 fm (0.4 fm). In reality, some of 
these small-size configurations are qq + gluons or qqqq, 
requiring a detailed model-dependent analysis. Even so, 
our result for the charge density places strong constraints 
on the pion's partonic structure at small distances. The 
study of dynamical models of PLCs in the pion and their 
comparison with the empirical charge density will be the 
subject of future work. 



V. LONG-RANGE PION STRUCTURE AND 
CHIRAL DYNAMICS 

To complete our study of the empirical transverse 
charge density in the pion we briefly want to comment 
on the possible role of chiral dynamics at large transverse 
distances. At b > m^ 1 = 1.5 fm the weighting factor 
KoiVib) in the dispersion integral Eq. ((3|) emphasizes the 
near-threshold region \ft — 1m v ~ few (see Fig. [TJd) , 
where the imaginary part of the form factor is governed 
by chiral dynamics and calculable from first principles In 
leading order of the chiral expansion, the imaginary part 
near threshold results from the pion loop graph with the 
7T7T 4-point coupling and is given by [38-40] 



7r _1 Imi^(t + «0) 



2 \3/2 



(t - Ami) 
6(47r/ w ) 2 Vi 



(12) 



Substitution of this result in Eq. ([3]) allows one to de- 
rive the leading exp(— 2m T 6) asymptotic behavior of the 
pion charge density at large distances; see Ref. [32| for 
details. Numerical analysis shows that the contribution 
from Eq. (jl~2j) to the charge density is negligible compared 
to the non-chiral density resulting from \ft ~ m p for all 
but the largest distances, reaching only ~ 30% of the dis- 
persion result at b = 2 fm. In the nucleon isovector form 
factor the chiral component of the charge density was 
found to become comparable to the non-chiral density 
at distances b ~ 1.7 fm (32j; the reason for its diminished 
role in the pion charge density is that the triangle grap h 
involving the ttN Yukawa coupling (see Fig. 1 of 32]), 
which gave the main contribution in the nucleon case, is 
absent for the pion. Account of higher-order chiral cor- 
rections does not substantially change the magnitude of 
the chiral component [loj]. We conclude that the trans- 
verse charge density in the pion is dominated by the p 
meson mass region for all distances of practical relevance, 
b < 2 fm. 



VI. SUMMARY AND DISCUSSION 

This paper shows how the pion form factor in the time- 
like region can be used to determine the transverse charge 
density. The timelike data greatly augment the mea- 
ger information available from spacelike pion form factor 
measurements, in particular in the region of high momen- 
tum transfers |i| > 1 GeV 2 conjugate to short transverse 
distances. Given the energy reach of the timelike form 
factor data, and the theoretical uncertainties involved 
in separating the real and imaginary parts, we estimate 
that p-n{b) is determined to an accuracy of ~ 10% at 
b = 0.1 fm, and substantially better at larger distances. 
The transverse density obtained from the full dispersion 
integral turns out to be surprisingly close to that ob- 
tained from a single zero-width p meson pole over a wide 
range. The empirical transverse density shows a strong 
rise at small distances, which points to a substantial pres- 
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ence of PLCs in the pion's partonic wave function and 
puts strong constraints on the pion GPD. 

In the work reported here we limited ourselves to a phe- 
nomenological analysis of the transverse density based 
on an existing parametrization of the timelike pion form 
factor data. Our results suggest several directions for 
further studies, both theoretical and empirical. 

The striking similarity of the empirical transverse den- 
sity to the simple vector meson dominance model over a 
wide region of b should have a dynamical explanation. 
Possible approaches to address this question are local 
quark-hadron duality or the dual resonance picture of 
QCD in the large- N c limit. 

The strong rise of the pion's transverse charge density 
at small distances calls for an explanation in terms of 
dynamical models of the pion's partonic structure. The 
key question is whether the required PLCs in the pion 
could be explained as the result of perturbative QCD 
interactions with large-size configurations, or whether 
non-perturbative interactions play an essential role. Of 
particular interest for addressing this question are mod- 
els which implement the non-perturbative short-distance 
scale associated with the spontaneous breaking of chiral 
symmetry in QCD such as the instanton vacuum model, 
which is known to give a reasonable description of the 
spacelike pion form factor at intermediate momentum 
transfers Q 2 ~ few GeV 2 [HI, S3- 

The dispersion result for the transverse charge den- 
sity at distances b ~ 0.1 fm depends sensitively on the 
phase of the pion form factor in the region of the lowest 
excited p states, Vt = 1 — 3 GeV. While the alternat- 
ing sign of the coupling of successive resonances is sug- 
gested by theoretical considerations, it would be worth- 
while to attempt independent experimental tests of this 
key assumption. This could be done through coherent 
photo- or electroproduction of two pions on nuclear tar- 
gets, which can be analyzed in the spirit of the gener- 



alized vector meson dominance model; see Ref. 43[ and 
references therein. Such measurements become feasible 
with the 12 GeV Upgrade of Jefferson Lab. 

The recent CLEO data [T^, which are difficult to ex- 
plain in the dual resonance framework commonly used to 
parametrize the high-energy region of the pion form fac- 
tor, may have a significant effect on the charge density at 
distances b < 0.1 fm. Confirmation of this experimental 
result and more data in the energy region \fi = 3 — 4 GeV 
would certainly be welcome. It would be interesting to 
explore ways to include these data in a dispersion analy- 
sis with more general parametrizations of the imaginary 
part. 

The new application of the timelike pion form factor 
described here once more underscores the importance of 
analyticity in relating observables measured in different 
kinematic regions. It would be helpful if phenomeno- 
logical parametrizations of the form factors such as (33j 
employed a framework which strictly respects analyticity, 
e.g. by using analytic functions like the GS form, or by 
parametrizing only the spectral strength on the physical 
cut and generating the real part by a dispersion integral. 
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Use of a subtracted dispersion relation in Eq. ([T| would 
lead to an expression for the transverse density which 
differs from Eq. <(3j by a delta function 5^ 2 \b). Subtrac- 
tions therefore have no influence on the dispersion result 
for the transverse density at finite b. 
Reference (3^| also provides an alternative fit in which 
the resonances are described by the BW form. As men- 
tioned above, the BW form has a spurious branch cut 
singularity at t = 0, which in the full GS form is can- 
celed by the off-shell term. This has the consequence that 
the imaginary part obtained with the BW parameters 
of Ref. [33( does not satisfy the normalization condition 
7r_1 O dt' lmF n (t')/t' = 1 exactly, but with a small 
discrepancy related to the spectral strength on the un- 
physical cut (the BW ansatz in [33| was normalized to 
unit value at t = 0, not to unit integral over the imag- 
inary part). In our dispersion analysis we therefore use 
the GS parametrization, which has the correct analyticity 
structure and is free of this problem. 
This argument assumes that the change of ImFV(t) does 
not substantially alter the cancellation of successive res- 
onance contributions implied by the dual resonance pic- 
ture. The effect of an "unbalanced" resonance will be 
estimated separately below. 

We note that the fit of Ref. [13] also does not reproduce 
the value of |fV| 2 at y/t = 3.1 GeV extracted from the 
J/ip — > tttt decay; see Ref. [33l ] for a critical discussion of 
this datum. 

The valence quark density in Eq. ()10|) is assumed to be 
normalized to unit integral, J dxq 7T (x) — 1. Note, how- 
ever, that Eq. (|10[) is intended as a model for laxge-x 
configurations and should be integrated only over this 
region. 



